Dielectric function of the semiconductor hole gas 
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The semiconductor hole gas can be viewed as the companion of the classic interacting electron gas 
with a more complicated band structure and plays a crucial role in the understanding of ferromagntic 
semiconductors. Here we study the dielectric function of an homogeneous hole gas in zinc-blende III- 
V bulk semiconductors within random phase approximation with the valence band being modeled 
by Luttinger's Hamiltonian in the spherical approximation. In the static limit we find a beating 
of Friedel oscillations between the two Fermi momenta for heavy and light holes, while at large 
frequencies dramatic corrections to the plasmon dispersion occur. 
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The interacting electron gas, combined with a homo- 
geneous neutralizing background, is one of the paradig- 
matic systems of many-body physics [ll-Q • Although ob- 
viously a grossly simplified model of a solid-state system, 
its predictions provide a good description of important 
properties of three-dimensional bulk metals and, in the 
regime of lower carrier densities, n-doped semiconductors 
where the electrons reside in the s-type conduction band. 

On the other hand, in a p-doped zinc-blende IITV 
semiconductor such as GaAs, the defect electrons or holes 
occupy the p-type valence band whose more complex 
band structure can be expected to significantly modify 
the electronic properties. Moreover, the most intensively 
studied ferromagnetic semiconductors such as Mn-doped 
GaAs are in fact p-doped with the holes playing the 
key role in the occurrence of carrier-mediated ferromag- 
netism among the localized Mn magnetic moments [J]. 
Thus, such p-doped bulk semiconductor systems lie at 
the very heart of the still growing field of spintronics @ , 
and therefore it appears highly desirable to gain a deeper 
understanding of their many-body physics. 

Following the above motivations, we investigate in the 
present letter the dielectric function of the homogeneous 
hole gas in p-doped zinc-blende IITV bulk semiconduc- 
tors within random phase approximation (RPA)[ll-Q. 
The single-particle band structure of the valence band 
is modeled by Luttinger's Hamiltonian in the spherical 
approximation In previous work we have studied 
the same system using Hartree-Fock (HF) approxima- 
tion fr]. A key result here is the observation that in 
a fully selfconsistent solution of the HF equations the 
Coulomb repulsion among holes modifies the Fermi mo- 
menta compared to the non-interacting situation. In par- 
ticular, the selfconsistent solution of the HF equations 
is not equivalent to first-order perturbation theory as 
it the case for the ordinary electron gas [ll-Q- More- 
over, we mention recent studies of the dielectric function 
in two-dimensional electron systems with spin-orbit cou- 
pling [i!, loj and two-dimensional hole systems 10]. Other 
recent related studies have dealt with the dielectric func- 
tion of planar graphene sheets where an effective spin is 



incorporated by the sublattice degree of freedom 11 , 13] . 

Luttinger's Hamiltonian describing heavy and light 
hole states around the T in III-V zinc-blende semicon- 
ductors reads (a 
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where niQ is the bare electron mass, p is the hole lattice 
momentum, and S are spin-3/2-operators. The dimen- 
sionless Luttinger parameters 71 and 72 describe the va- 
lence band of the specific material within the so-called 
spherical approximation The above Hamiltonian is ro- 
tationally invariant and commutes with the helicity op- 
erator A = {k ■ S)/k, where k — p/h is the hole wave 
vector. Thus, the eigenstates of ([Ij can be chosen to be 
eigenstates of the helicity operator with the heavy (light) 
holes corresponding to A = ±3/2 (A = ±1/2). The en- 
ergy dispersions are given by Sh/i{k) = h^k"^ /2mh/i where 
Tnh/i = Wo/ (71 T 272) is the effective mass of heavy and 
fight holes, respectively. 

Combining the above single-particle Hamiltonian with 
Coulomb repulsion among holes and a neutralizing back- 
ground, the dielectric function within RPA is generally 
given by 
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where V{k) is the Fourier transform of the interaction 
potential, and the free polarizability reads 
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Here /(fc. A) are Fermi functions, and the explicit form of 
the four-component eigenspinors IyaI^)) of the Hamilto- 
nian ([l} has been given in Ref. [7|]. The mutual overlap 
of these eigenspinors entering the above expression is a 
key feature of the semiconductor hole gas. 
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In general, an exact evaluation of the free polarizability 
([3]) is, even in the limit of zero temperature, a formidable 
task and clearly more complicated than the case of the 
spinless electron gas. Therefore we shall be content here 
with zero-temperature properties concentrating on the 
static limit, and on the regime of large frequency and 
small wave vector. In the former case (w = 0) an already 
quite tedious calculation yields 
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where kh/i = ^j2mhiiSF 1^ are the Fermi wave numbers 



for heavy and light holes at Fermi energy ef- 
called Lindhard correction L is given by 
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and the function H is defined as 
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Remarkably, one can express the polarizability entirely in 
terms of the arguments k/2kh, k/2ki, and k/kh + ki with 
the latter one being the harmonic mean of the two former. 
In the limit rUh = mi (i.e. kh — ki =: kp) one obtains 
the usual result Xo(k,Q) — —D{ep)L{k/2kp) for charge 
carriers without spin-orbit coupling where D{e) is the 
density of states [13| . The full polarization ^ at nih ^ 
mi, however, has a clearly mor complictcd structure. 

On the other hand, considering Coulomb repulsion, 
Vik) = / SrS^k"^ , and using the long- wave approxi- 
mation xo{k,Q) ~ xo(0,0) leads to the usual Thomas- 
Fermi (TF) screening, 0) w 1 - k'^p/k'^ with 
kxp = (e^/ereo)3n/2£_F. Here Sr is the background di- 
electric constant taking into account screening by deeper 
bands, and the hole density is given by n — Uh + ni, 
nh/i 
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The full screened potential of a pointlike probe charge 
Q is given by 
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whose asymptotic behavior is determined by the singu- 
larities of the integrand and its derivatives IJ] . Here the 
first derivative has singularities at A: = 2kh and k = 2ki 
while at A: = kh + ki all singular contributions cancel out. 
As a result, the Lighthill theorem [31 yields for large 
distances r 
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where 



cos(qr) 



(9) 



and oo — 'ineoh? /mQe'^ being the usual Bohr radius. 
Thus, we observe a beating of Friedel oscillations be- 



The so- tween the two wave numbers 2kh/i- Note that, differ 
ently form the expression for the dielectric function it- 
self, the wave number k = kh + ki does not occur in 
the Friedel oscillations since the non-interacting ground 
state of the hole gas has singularities in the occupation 
numbers at A: = kh/i but not at k — {kh + fc;)/2. Fig. [T] 
shows the Friedel oscillations according to Eq. ^ along 
with a numerical evaluation of the full Fourier integral ([7]) 
for p-doped GaAs with a hole density of n = 10^'^cm^'^, 
which is a very typical value for Mn-doped GaAs [4] . One 
might argue whether one should replace the Fermi mo- 
menta kh/i with renormalized values arising from a fully 
self-consistent solution to the HF equations. However, 
at large densities this renormalization becomes negligi- 
ble 0. 

The beating of Friedel oscillations illustrated in the 
figure is a peculiarity of the holes residing in the p-type 
valence band and should be observable via similar scan- 
ning tunneling microscopy techniques as used in metals 
fisj and n-doped semiconductors [l^. Moreover, as theo- 
retical studies have revealed, such oscillations can have a 
profound impact on the magnetic properties of ferromag- 



netic semiconductors III, ll8|. Moreover, Fig. [T] shows 



the amazing accuracy of the asymptotic expression ([8]) 
obtained from the Lighthill theorem. 

Let us now turn to the regime of large frequencies and 
small wave vectors. Following Ref. [2\ we expand the 
denominators in Eq. ([3]) assuming huj >> eh/i{k) and 
hjj » {hkh/i/mh/i)hk. Within the two leading orders 
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FIG. 1: Friedel oscillations resulting from a numerical evalua- 
tion of the Fourier integral ([7| , and estimated via the Lighthill 
theorem (cf. Eq. (|8])) for p-doped GaAs with a hole density of 
n — lO^^cm"^. The inset shows the data at larger distances 
on a smaller scale. 
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For mil = mi the first three lines of the above expression 
reproduce again the standard textbook result j2| while 
all other terms vanish in this limit. On the other hand, 
if m/i ^ mi , one has contributions in order 1 /w'* that are 
independent of the wave vector k. Such terms are absent 
in the case of the standard electron gas where the contri- 
butions of order l/w^" are at least of order in the 
wave vector Q • The technical reason why such contribu- 
tions are present for the hole gas is that the expression 
ex2ik' + k)- eAi(fc') in Eq.© contains for |Ai| ^ IA2I 
an additive term which is independent of k (and vanishes 
for TO/,, = mi). These prima vista unexpected contribu- 
tions to the high-frequency expansion of the dielectric 



function will also occur in even higher orders. However, 
even in the two leading orders given in Eq. PH)) . they 
strongly modify the plasmon dispersion determined by 
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where the zero-order plasma frequency is given bv[l9| 
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and the dimensionless and density-independent coeffi- 
cients zi, V, w are given by 
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Clearly, the coefficients u and w vanish for to,, = to, while 
from V one recovers usual textbook result for an electron 
gas without spin-orbit coupling Q. By expanding the 
square root in Eq. (jlip we have neglected higher contri- 
butions both in wave vector and in the density parameter 
n^/'^flQ cx {e F / hid^p'^ which is consistent with consider- 
ing only the first two leading orders in Eq. ([TU| . In fact. 
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AlAs 


0.47 


0.18 


10.0 


0.38 


17.7 


21.5 


-16.3 


AlSb 


0.36 


0.13 


12.0 


0.36 


49.7 


37.1 


-29.5 


GaAs 


0.5 


0.08 


12.8 


0.16 


195.4 


99.4 


-100.5 


InAs 


0.5 


0.026 


14.5 


0.052 


861.4 


451.9 


-473.1 


InSb 


0.2 


0.015 


18.0 


0.075 


1796.9 


919.2 


-958.8 



TABLE I: Material parameter and coefficients u, v, w of the 
plasmon dispersion (|12p for various III-V semiconductors. 

for usual p-doped bulk semiconductors n^^^ao is small, 
and to consistently obtain contributions to the plasmon 
dispersion being of higher order in the density would re- 
quire to extend the expansion (I10[) also to higher orders, 
which is computationally increasingly tedious and will 
lead to even lengthier expressions. Note that the disper- 
sion coefficients u, v, w depend entirely on material pa- 
rameters. In tableUwe have listed their numerical values 
for several prominent III-V semiconductor systems. As 
seen there, the coefficient u is remarkably large leading to 
a substantial enhancement of the long- wavelength plasma 
frequency w^(0) — {lu^^)'^{1 + u{n^^^ao), even at small 
densities, compared to the naive guess w^(0) « (wp""*)^. 
On the other hand, v and w differ in sign and are of 
quite similar magnitude resulting in a dramatic flatten- 
ing of the plasma dispersion compared to the standard 
case rrih = mi where w vanishes. Moreover, the sum 
V + w can even become negative leading to a plasmon 
dispersion bending downwards around zero wave vector. 
In fact the sign of v+w is entirely determined by the ratio 
mi/mfi where negative values occur for mi/mh ^ 0.18. 
Remarkably, GaAs lies very close this threshold showing 
already such a qualitative change in the plasmon disper- 
sion. This trend is further enhanced in the cases of InAs 
and InSb. 

In summary, we have studied the dielectric function of 
the homogeneous hole gas in p-doped zinc-blende III-V 
semiconductors. In the static limit we predict additional 
beatings of the Friedel oscillations which should be ex- 
perimentally detectable via state-of-the-art scanning tun- 
neling microscopy. At high frequencies and small wave 
vectors the plasmon dispersion gets dramatically altered 
compared to the textbook case of the usual electron gas. 

I thank J. Repp for useful discussions and Deutsche 
Forschungsgcmcinschaft for support via SFB 689. 
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